Optimal path for a quantum teleportation protocol in entangled networks 
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Bellman's optimality principle has been of enormous importance in the development of whole branches of ap- 
plied mathematics, computer science, optimal control theory, economics, decision making, and classical physics. 
Examples are numerous: dynamic programming, Markov chains, stochastic dynamics, calculus of variations, 
and the brachistochrone problem. Here we show that Bellman's optimality principle is violated in a teleportation 
problem on a quantum network. This implies that finding the optimal fidelity route for teleporting a quantum 
state between two distant nodes on a quantum network with bi-partite entanglement will be a tough problem and 
will require further investigation. 



Finding the route between two nodes of a given graph such 
that the sum of the weights associated with the links within 
this path is minimized is arguably one of the most fundamen- 
tal problems in graph theory (]]]. Decades ago, an algorithm 
to solve it was proposed by Dijkstra [2]. Its applicability re- 
lies on the ability to compute the optimal path for a large-scale 
network based on the optimization performed at smaller parts 
of it, with initial conditions given by other parts. As stated by 
Dijkstra, to find the path of minimal total length between two 
given nodes P and Q, "we use the fact that, if R is a node on 
the minimal path from P to Q, knowledge of the latter implies 
the knowledge of the optimal path from P to R". This can be 
understood as a particular instance of the more general opti- 
mality principle by Bellman [3[[: for problems satisfying this 
principle, the global optimal solution can be determined in 
terms of local optimal ones for smaller subproblems (optimal 
substructure). This property has been extremely useful in the 
study of dynamic programming, control theory, economics, 
and Markovian stochastic processes. 

Several problems originally thought in classical scenarios 
have been subsequently redefined and studied when dealing 
with quantum systems. Significant differences in the pre- 
dictions are found when focusing on atomic and sub-atomic 
scales, but even some macroscopic phenomena can be under- 
stood only through a complete analysis by means of quantum 
theory. In particular, the role of entanglement [4] in different 
processes in nature is currently under deep investigation |5|]. 
On the other hand, the generalization of classical information 
theory to quantum scenarios has paved the way to the develop- 
ment of quantum information theory, with rules that are fun- 
damentally different from classical ones @]. 

Here we focus our attention on a particular scheme of quan- 
tum information theory, known as quantum teleportation, that 
has proved to be one of the most striking applications of 
entanglement as a resource |7|]. By means of a shared en- 
tangled channel and a conditional local operation, a receiver 
(Bob) is able to reconstruct the unknown state of a qubit given 
to a sender (Alice), after she performs a joint measurement 
and communicates classically her result to him. Consider- 
able effort has been made in demonstrating quantum telepor- 
tation experimentally by means of polarized entangled pho- 
tons HI], squeezed-state entanglement [9], liquid-state nuclear 



magnetic resonance ifioll . and trapped ions ifTHl . In this paper, 
we show that Bellman 's optimality principle is violated when 
we consider the teleportation protocol on a quantum network. 

Although further generalizations and extensions of Dijk- 
stra' s algorithm have been studied [12], here we will shortly 
illustrate its working principles. If we want to find the shortest 
path (in terms of the sum of the weights) between two nodes 
A and B of a graph, we need first of all to assign a number 
to each node. We give the value to A and +oo to all the 
others. Then we start from A (which we define as the "incom- 
ing" node, until it is marked) and we consider its neighboring 
nodes. To each of them, we assign the minimum between its 
current value (in this case +oo) and the sum of the value of 
the incoming node (in this case 0) and the weight of the link 
between the latter and the node under investigation. After we 
have done this for all the neighboring nodes, we mark A, and 
we no longer need to consider it. We then focus our atten- 
tion on the unmarked node in the network with the smallest 
number, and we highlight the link that "contributed" to gen- 
erate its actual value. We proceed in the same way with this 
new unmarked node taking the role of A in the previous de- 
scription. The process stops when there is a highlighted link 
ending on B. The only highlighted path going from A to B is 
the shortest path between these two nodes of the graph. 

The investigation on quantum teleportation has been mainly 
focused on the case where Alice and Bob directly share an 
entangled resource. Interesting results have been obtained 
when a sequence of teleportations has to be performed along 
a chain of nodes 111311 . In order to provide a broader scenario 
for quantum teleportation, we consider here a more complex 
setting: Alice and Bob are connected through a network in 
which each node shares entangled channels (that are, in gen- 
eral, non-maximally entangled) with others. A sketch of the 
situation under investigation is presented in Fig.QJa)- This is 
in line with the recent interest that arose in the scientific com- 
munity about the possibility to realize quantum networks ffl . 
Networking distant nodes is a fundamental step in design- 
ing and building distributed quantum computers, as well as 
in the implementation of large-scale highly secure quantum 
communication protocols. This concept paved the way for the 
promising idea of quantum internet [15]. We focus, in partic- 
ular, on quantum networks with just bi-partite entanglement 
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FIG. 1 : (a) Sketch of the general setting considered: Alice and Bob 
are connected through a network in which each node shares entan- 
gled channels with others, (b) Simplified scheme: Alice can teleport 
her state to Bob's position either by means of the direct link A — B 
(path Vi) or through the link A - C and then C — B (path V2). 



(therefore, there is no multi-partite entanglement in our net- 
work). 

The scenario described here resembles the shortest-path 
problem aforementioned. Surprisingly, we will find that Dijk- 
stra's algorithm can be adapted for solving our problem, while 
the entangled channels are pure. When the nodes in our net- 
work are linked by means of mixed resources, we can still 
consider the analogy with the graph theory problem, but in 
general we can no longer use Dijkstra's method. Moreover, 
we will find that our problem does not satisfy the broader Bell- 
man's principle and then also all the other algorithms based 
on it cannot be used. We will discuss how to determine for 
which cases the principle is still satisfied. 

Let us start considering a simple scenario, where Alice, 
Bob, and Charlie (an agent located in a third node of the net- 
work) are linked by means of non-maximally entangled pure 
states, as shown in Fig. [Tib). The extension to more complex 
network structures and mixed channels is discussed later. We 
want to stress again that we are dealing with the case of bi- 
partite entangled resources here. This assumption, together 
with the others within this paper, is a case that can be exper- 
imentally tested with current state-of-the-art quantum tech- 
nologies. A general non-maximally entangled pure state of 
two qubits can be cast, in an appropriate basis, as |$( e ))i 2 = 
cos 6> |0)i |0) 2 + sin0 1) i|l) 2 , where 9 e [0, w/4] is a param- 
eter depending on the amount of entanglement present in the 
channel (related to Schmidt coefficients) 111 611 . As a measure 
of entanglement, we will use the negativity 11711 of the density 
matrix p\,2 describing the state of qubits 1 and 2 (in the case of 
pure states, this corresponds to pi.2 = |$^)i,2(^ ^ |). In the 
scenario considered so far, we have that the negativity A/i.,2 of 

the state is simply A/1,2 = sin 29. 

Alice has two possible choices to teleport her input state to 
Bob's position (we assume that the standard teleportation pro- 
tocol is used at each step): either she can use the direct link 
A — B [path V\ in Fig. QIb)], or she can teleport her input 
state to Charlie through the link A — C, and he will subse- 
quently teleport it to Bob through the link C — B [path V2 
in Fig.QJb)]. Let us define the negativity of the three shared 



entangled states, corresponding to these three links, as Na,b, 
Na,c an d Mc,b, respectively. The result of a deterministic 
teleportation protocol (we are not dealing with probabilistic 
schemes lfl8ll in this paper) is the final state p^f at Bob's 
location. This is, in general, a mixed state that is not exactly 
the same as the initial one = \^)a{4'\ mat Alice wanted 
to teleport. Although we consider here a pure input state, the 
extension to the mixed case is straightforward. Moreover, we 
focus the investigation first on azimuthal states of the form 



a\0) + /3|1), with \a\ 2 + |/3| 2 = 1 and a, /3 both real llHI^. 

In order to evaluate the performance of the scheme, we need 
to estimate how close pg is to /0;J by means of the tele- 
portation fidelity T = a{iP\Pb W)a- This value clearly 
depends on the state to transfer; therefore an average over 
all the possible input states is in order. Let us define T as 
the average of T over all the possible \iP)a uniformly dis- 
tributed. Interestingly, the average fidelity if Alice performs 
a direct teleportation through the link A — B can be simply 
cast as Ta-b = (3 + A(a,s)/4. On the other hand, if she 
chooses to use the two-link path A — C — B, the fidelity is 
Ta-c-b = (3+ Ma, cA/c,b)/4 (a similar result has been ob- 
tained in Ref. J2JJ]). Clearly, the general case through a more 
complex path V can be written as Tp = (3 + Yiiev^)/^ 
where A/j denotes the negativity of the entangled channel cor- 
responding to the link i within the path V . Here it is easy to 
note how to map the maximization of the teleportation fidelity 
to the aforementioned shortest-path problem. If we consider 
the value - In A/j (0 < A/i < 1, so we have -InA/j > 0) 
as the weight of link i, the maximization of Tp corresponds 
to the shortest path solution (i.e., to finding the path V that 
minimizes — J^iev In A/j). 

Let us now extend the analysis to the case where the shared 
resources are not pure. We first start by describing the main 
result, obtained for the interesting class of AT-states B22I1 . This 
includes, among others, maximally entangled Bell states and 
Werner states B23I1 . This class is relevant in many physical 
settings, and, for this reason, it has recently attracted the at- 
tention of the scientific community and has been investigated 
theoretically as well as experimentally [24]. An AT-state of 
two qubits is described by the density matrix 



fan aiA 

a 2 2 a 2 3 

a 32 a 33 

\a^i a^J 



(1) 



with ajk = a^j (ajy denotes the complex conjugate of akj). 
Following our discussion, we study the average fidelity when 
Alice teleports her input state through a general path. Straight- 
forward but cumbersome calculations lead to 
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Here a-^ is the corresponding element of the density matrix 
describing the entangled state of link i within the path V. It 
is now straightforward to see why, in general, Dijkstra's algo- 
rithm can no longer be used to solve the problem of finding 
the corresponding shortest path: the presence of two different 
non-constant terms in the formula for T-p in Eq. (f2|, each one 
equal to the product of values associated with the links in the 
path, does not allow us to assign to the links a single weight 
with the additive property. We can show more neatly the vio- 
lation of Bellman 's optimality principle with the following ex- 
ample: assume that, by applying the fidelity formula in Eq. (O 
for the case of X-states, we determine that V\ is the optimal 
path to go from A to B in Fig. 02b). Next, imagine that we 
add a node D connected only to B through an entangled link 
of the same class. If we redefine our task as teleporting the 
unknown state from A to D, then it is impossible to guaran- 
tee that the optimal path for this new problem will contain the 
subpath V\ (but, according to Dijkstra's protocol, if the opti- 
mal path from A to D contains V2 and not V\, then V\ cannot 
be optimal from A to B). In other words, the optimization 
depends not only on the present node and the total "distance " 
accumulated before but also on the whole structure of the net- 
work. This justifies why finding the optimization in this type 
of networks is not possible by means of any algorithm based 
on Bellman's optimality principle. 

Nevertheless, the problem of maximizing the teleporta- 
tion fidelity can still be mapped into finding the shortest 
path in a graph, but this time each link i has two different 
weights and vi and the distance along a path is obtained 
as riiGP + Ilie-p v i + Vt where r\ is a constant. The 
corresponding weights are ^ = a[ 1 ^ — a^l — 033 + a^l 
and Vi = a^l + + a 3 2 + a^i, respectively. When 

a 22 = a 33 = 0' we nave a n + a 44 = 1 (due t0 tne f act 
that the trace of a density matrix is always equal to 1), and 
Dijkstra's algorithm can be used again. This condition means 
that the shared entangled state has to be confined in the sub- 
space spanned by {|00), |11)}. 

We want to stress here that clearly this result is not a pure 
quantum effect but, as stated above, it depends on the fact that 
the teleportation fidelity in the investigated scenario cannot be 
related to individual link weights with a pure additive or mul- 
tiplicative property. This can be, of course, also the case in 
particular problems on classical networks, in which the figure 
of merit that we want to maximize (or minimize) behaves nei- 
ther additively nor multiplicatively under the composition of 
two channels within the network: Bellman's optimality prin- 
ciple is violated, and all the algorithms based on it cannot be 
used. However, it is interesting that the scenario presented 
here gives an explicit and intuitive example of this fact in a 
quantum case. 

Let us now consider a significant example of an X-state, 
a model that can be easily seen as a generalization of two- 
qubit Werner states [23J, in which \$( s ') takes the place of 
the Bell pair 112011 . In this case, the non-maximally entan- 
gled mixed state of qubits 1 and 2 is represented by the den- 



sity matrix p lf2 = p w \$ {e) ){® {e) + (1 - Pw) I (4) / 4 , with 
Pw S [0, 1] and I" being the 4x4 identity matrix. This can 
be, for instance, the result of an initially pure general entan- 
gled state (represented by that is successively affected 
by isotropic random noise due to the interaction with the en- 
vironment |0. Let us define the state of qubits j and k as 
p. ife =p^ fe )|$(^))($(%,,)| + (l-^ fe ))l(4)/4, with^fc 
being the angle in the definition of |$( fl 3. fe )}. The average tele- 
portation fidelity in this setting is 

?v = \[2 + + II Pw si*(20;)]. (3) 

This shows that, for this model of mixed resources, Dijkstra's 
algorithm (and, more generally, Bellman's optimality princi- 
ple) can no longer be used to solve the problem of finding the 
shortest path corresponding to the maximization of the aver- 
age teleportation fidelity. 

Now, let us comment on the restriction we have considered 
on the initial state to be teleported. Clearly, its more general 
form a\0) + is without any condition on a and (3, apart 
from \a\ 2 + \f3\ 2 = 1 (they are complex numbers). The set 
of all the possible initial states over which we have to take the 
average is thus different. However, if we average over this set, 
the formulas for the fidelity are cumbersome and cannot be 
cast in a simple form as for the case investigated above. For 
the sake of completeness, a longer technical work, including 
the formulas for a general initial state and further details, is 
in preparation. However, we want to anticipate here that, also 
in that case, Bellman's optimality principle is violated for the 
same conditions as those discussed in this paper. 

The formulas presented so far allow one to obtain the max- 
imum teleportation fidelity when the strategy to follow con- 
sists only in deciding an optimal route and performing single 
teleportation steps for each link of the path. However, one 
can outperform this value of fidelity including a probabilis- 
tic preparation stage in the scheme. In contrast to probabilis- 
tic teleportation [18], the teleportation proposed here is still 
deterministic. The probabilistic nature is present only in the 
preparation stage, in which the network is "arranged" in a way 
so as to achieve the best fidelity. After the preparation stage, 
one can perform the sequence of teleportations and obtain de- 
terministically Bob's final state. The fidelity will depend on 
the result of the preparation stage and its minimum value will 
correspond to the maximum fidelity calculated on the original 
graph according to the formulas presented in this paper. 

Let us explain the method by considering again the easy 
instance of network sketched in Fig. Htb). Suppose that the 
shared resources are pure and Na-b > Na-cNc-b, with 
the generalization to more complex network structures and 
mixed channels being straightforward. According to the pre- 
vious discussion, the optimal path should be the direct link 
A — B. The entangled resources corresponding to the links 
A — C and C — B are thus not needed for the implemen- 
tation of the teleportation protocol in order to achieve the 
fidelity Ta-b = (3 + A/a-b)/4. We can therefore use 
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the couple of entangled pairs A — C and C — B to per- 
form a probabilistic entanglement swapping and concentra- 
tion scheme in order to create an additional link A — B. For 
instance, by measuring the two qubits at location C on a 
Bell basis, one can obtain a new entangled resource A — B 
with negativity N' A _ B = 2Ma-cNc-bIi and a success 
probability p = (1 - ^Jl - Wa-c\/1 = A/c-s)/2. Here 
7 = [2— cos(#i— 5a)— cos(5i+£a)]. 5i = (arcsiri7v0t-C')/2, 
and (52 = (arcsinA/c 1 — s)/2 (clearly, this is only an exam- 
ple of the swapping-concentration schemes that can be per- 
formed: by properly choosing the basis for the measurement 
at location C, it is possible to have a trade-off between the fi- 
nal negativity of the new link A — B and the success probabil- 
ity). If the swapping-concentration scheme has been success- 
ful, the use of this link provides a higher teleportation fidelity 
than the one that could be obtained by means of a double tele- 
portation through the two original links A — C and C — B 
and can also outperform the one through the original link 
A — B. In this way, the average fidelity of teleportation has 
been increased. In the case where the swapping-concentration 
scheme has not been successful, both the state to teleport and 
the original link A — B are unaffected, therefore one can still 
use the original path in order to achieve the maximum fidelity 
according to the formulas without the preparation stage. 

The simple quantum scenario described here poses a very 
interesting problem. Once large quantum networks based on 
bi-partite entanglement are built up, their optimization during 
regular exploitation will require the development of new al- 
gorithms. Even though the optimal formulas we have derived 
can be beaten by allowing the protocol to become probabilis- 
tic, the violation of Bellman's optimality principle cannot be 
circumvented in the case where the probabilistic entanglement 
swapping for the nodes not belonging to the optimal route 
fails. 
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